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Abstract We use the redshift Hubble parameter H{z) data derived from relative galaxy 
ages, distant type la supernovae (SNe la), the Baryonic Acoustic Oscillation (BAO) peak, 
and the Cosmic Microwave Background (CMB) shift parameter data, to constrain cos- 
mological parameters in the Undulant Universe. We marginalize the likelihood functions 
over h by integrating the probability density P oc 1"^ . By using the Markov Chain 
Monte Carlo (MCMC) technique, we obtain the best fitting results and give the confi- 
dence regions on the 6 — ilmo plane. Then we compare their constraints. Our results show 
that the (z) data play a similar role with the SNe la data in cosmological study. By pre- 
senting the independent and joint constraints, we find that the BAO and CMB data play 
very important roles in breaking the degeneracy compared with the ll{z) and SNe la data 
alone. Combined with the BAO or CMB data, one can improve the constraints remark- 
ably. The SNe la data sets constrain fimo niuch tighter than the Hiz) data sets, but the 
}i{z) data sets constrain h much tighter than the SNe la data sets. All these results show 
that the Undulant Universe approaches the ACDM model. We expect more ll[z) data to 
constrain cosmological parameters in future. 
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1 INTRODUCTION 



Many cosmolo g ical ob servations jBennett et al.l (l2003l); Spergel et al.l (l2007h; Eisenstein et al. [i2005h . 
iKowalski et"ai1 (|2008|) ). such as distant SNe la jRiess et all (Il998l) ; iPerlmutter et al.1 (Il999l)). the 



subseq uent CMB measurement by Wilkinson Microwave Anisotropy Probe (WMAP) dSpergel et al 



2003[)) and the large scale str ucture survey by Sloan Digital Sky Survey (SDSS) dTegmark et al 



( l2004al) ; iTegmark et alj ( l2004bh ) etc. support that our Universe is undergoing an accelerated expan- 
sion, and these data hav e mapped the universe s ucce s sfully b y constrain i ng cosmological parameter s 
(IPerlmutter et alj ( Il999h; iTegmark et al] (|2000|); Izhul (|2004|); Izhan et al.1 ( |2006|) ; iLineweaveil ( Il998h ; 
lEfstathioul dfgggh ; IZhan et al.l (l2005h : iMaartens et al.l d2006t) '). Meanwhile, many different observa- 
tions, such as the redshift-dependent quantities, are used usually. The l uminosity distance to a par- 
ticular class of objects, e.g. SNe la and Gamma-Ray Bursts (GRBs) dNesseris & PerivolaropoulosI 
d2004l) ). the size of the BAO peak dete c ted in the large-scale correlation function of luminous 
red galaxies from SDSS dEisenstein et al.l d2005l) ) and the CMB data obtained from the three-year 
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WMAP estimate dWang & Mukherjed (l2006h : ISpergel et all (l2007l) ) are redshift-dependent. The BAO 
and CMB data have been combined to constrain the cosmological parameters widely. Recently, 
the Hubble expansion rate at different redshifts can be obtained from the measurement of relative 
galaxy ages directly. A s a fun c tion of redshift z, t he H(z) data have been u s ed to test cosmologi - 
cal models dGong et alj (l2008l): lYi & Zhand (|2007|): ISamushia & Ratral (l2006h: IWei & Zhan'S (l2007h: 
IZhang et al.' (lOOf): 'Wei & Zhang' fe OOTb: IZhang & Wul (|2007|): iDantas et alj ( l2007h : IZhang &Zhul 
T2OO7); Wu & Yu (2007 ); Wei & ZhanJ "( l2007h ). Using different data to constrain parameters can pro- 
vide the consistency c hecks be t ween e ach other, and the combination of d ifferent data can also make 
the constraints tighter dLi et alJ (Eooi i lZhao et al.l (l2007h : IXia etaP (l2006h l. So far, few work has been 
done in comparison. We'd like to study the comparison of the SNe la, H jz), BAO and CMB d ata in 
this letter A lthough lots of cosm ological models, e.g. the Quintessence ( Caldwell et al] (Il998h ). the 
brane world dOeffavet et al. (|2002|) ). the Chaplygin Gas dAlcaniz et alj d2003l) ) and the holographic dark 
energy models dKe & Lil d2005h ). are extensively explored to explain the accelerati on of the univers e. 
ACDM model, which is a standard and popular model, has been studied by Lin et al. dLin et al] d2008l) ). 
In fact, there are some problems in the standard model, and this model is not consistent with the real 
universe. Scientists propose d a Undu l ant Un iverse, which is characterized by alternating peri ods of ac- 
celeration and deceleration ('Gabriela'('2005')). This model can remove the fine tuning proble m dVilenkinI 
d200lh : IJaume et al.l Jl999) ; Bludman aOOCj); I Garriga & VilenkinI d200ll) : [StewartI J2OOOI) ). and solve 
the coinci dence scandal between t he observed vacuum energy and the current matter density, with some 
details in dPeebles & Ratral d2003b ). The equation of state (EOS) of the Undulant Universe is 



a; (a) = — cos(blna), 



(1) 



where the dimensionless parameter h co ntrols the frequen c y of t he accelerating epochs. Due to most 
inflation models predicting fik < 
Hubble parameter is 



10-5 dTegmark & ReesI dl998h ). we assume spatial flatness, so the 



H'^{a) ^ Hl{Q.n^aa + f^Aoa ^ cxp[- sin(61na)]), 



(2) 



where f^mo, f^Ao are the density parameter of matter and dark energy, respectively. Furthermore, 
S. Nesseris and L. Perivolaropoulos have enhanced the fact that an oscillating expansion rate 
ansatz ha s provided the best fit to the data among many ansatz by using recent supernova 
data CNesseris & PerivolaropoulosI (|2004|) ). 

In this letter, we concentrate on the comparison of different data sets on the Undulant Universe. In 
§11, we simply introduce these data sets that we use. Then the constraints are shown in §111. Finally, we 
show some important conclusions and discussions in §IV. 



2 OBSERVATIONAL DATA 
2.1 SNe la data 

There has been some calibrated SNe la data with high confidence in (iKowalski et al.l d2008l) ). This 
data set contains 307 data and the redshifts of these data span from about 0.01 to 1.75. These data give 
luminosity distances c^l [zi] and the redshifts z, of the corresponding SNe la. In a Friedmann-Robertson- 
Walker (FRW) cosmology, considering a flat universe, the luminosity distance is, 

dL(z) = £il±ilj-(z), (3) 
^0 

The function J-{z) is defined as J-{z) ~ dz/E{z), with E{z) — H{z)/Hq = 
\JriMoa~^ + r^Aoa-^ 6xp[| sin(&lna)]. Hq is the Hubble constant. The distance modulus is 

^(z) = m - M = 5 log = 42.39 + 5 log ^^-^(z), (4) 
lOpc n 
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with h = i/o/lOO km s ^ Mpc ^, where ni and M are the apparent and absolute magnitudes, 
respectively. 



2.2 The BAO Data 



As the acoustic oscillations in the relativistic plasma of the early universe will also be imprinted onto the 
late-time power spectrum of the non-relativistic matter (Eisenstein & Hu ( 1998)), the acoustic signatur es 
in the large-sca le clustering of galaxie s yield additional tests for cosmology dSpergel et alj (|2003|) ). 
Eisenstein et al. jEisenstein et alJ (l2005h ) choose a spectroscopic sample of 46748 luminous red galaxies 
from SDSS to find the peaks successfully. These galaxies cover 3816 square degrees and have been 
observed at redshifts up to z = 0.47. The peaks are described by the model-independent ^-parameter, 
yl-parameter is independent of Hubble constant Hq, 



-F2(2l))]l/^ (5) 



where zi = 0.35 is the redshift at which the acoustic scale has been measured. Eisenstein et al. suggested 
that the measured value of the ^-parameter is ^ = 0.469 ±0.017 jEisenstein et al.l (l2005l) ). 




2.3 The CMB Data 



The expansion of the Universe has transformed the black body radiation, which is left over from the B ig 
Bang, into the nearly isotropic 2.73 K cosmic microwave background (CMB) dBernardis et al.l ( |2000|) ). 
The whole shift of the CMB angular power spectrum is determined by the CMB shift parameter TZ. The 
shift parameter perhaps the most model-independent parameter , is also independent of Hq. It can be 
derived from CMB data ( iBond et alJ ( ll997h : IOdman et alj ( l2003h ). 

n = v^-^(2r), (6) 

where = 1089 is the redshift of recombination. Using the Markov Chain M onte Carlo (MCMC) 
chains from the analysis of the three-year results of WMAP jSpergel et al.l (l2007l) ). Wang & Mukherjee 
compute the CMB shift parameter to get TZ = 1.70 ± 0.03 and demonstrate that its measured value is 
mostly independent of assumptions about dark energy (IWang & Mukherjeel(l2006h ). 



2.4 the H{z) data from relative galaxy ages 



The Hubble parameter H(z) data can be derived from the derivative of redshift 
cosmic time t, i.e. dz/dt Uimenez & Loeb (2002)), 



Hiz) = 



1 dz 

1 + z dt ■ 



with respect to the 



(7) 



Therefore, an application of the differential age method to old elliptical galaxies in the local universe can 
determine the value of the current Hubble constant. This is a direct measurement for H{z) through a 
determination of dz/dt. Jimenez et al. (Jimenez et al.l (l2003b ) have applied the method to a z < 0.2 
sample. Paying careful attention to uncertainties in the distance, systematics, and model uncertain- 
ties, they have selected Monte-Carlo techniques to demonstrate the feasibility by evaluating the er- 
rors. They also use many other different methods and yield consistent value for the Hubble constant. 
Hence, the H{z) data are reliable. With the availability of new galaxy surveys, it becomes possible 
to determine H{z) at z > 0. By using the differential ages of p assivel y evolving galaxies determined 
from the Gemini Deep Deep Survey ( 'GDDS ) ([A braham e t alj (20041) ) and archival data (iTreu et al.l 



J ) dXi - ■ 

(l200lb : iTreu et al.l (|2002|) : |Nolan et al.' (i2003al):lNolan et al.l (l2003bl)). Simon et a l deriv e a set of g( , 
data, w hich is listed in TablefTl cSimon et al.l ( 2005 ). see also Samushia & Ratra (l2006l) : iJimenez et al.l 
(120031) ). The estimation method in detail can be found in the work dSimon et al.l (120051) ). As z has a 
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Table 1 the iJ (z) data in units of km s ^Mpc ^ (ISimon et al (see 
ISamushia & Ratral (|2006|) : iJimenez et all (|2003|) also) 



z(redshift) 
H{z) 
la error 



0.09 
69 
12.0 



0.17 
83 
8.3 



0.27 

70 

14.0 



0.40 
87 
17.4 



0.88 
117 
23.4 



1.30 
168 
13.4 



1.43 
177 
14.2 



1.53 
140 
14.0 



1.75 
202 
40.4 



4 - 



2 - 



1 1 1 1 1 1 1 1 


H(z) only 








\ \ 








\ 

\ \ 




\ 

\ \ 
1 



0.2 



0.4 



0.6 




Fig. 1 The contour maps of 6 — i7niO for different data sets. The left panel is H{z) only, and 
the right panel is SNe la only. Confidence regions are at 68.3%, 95.4% and 99.7% levels from 
inner to outer respectively, for a flat Undulant Universe without a prior of h. 



relatively wide range, 0.1 < z < 1.8, these data are expected to provide a full-scale description for the 
dynamical evolution of our universe. The application of the H(z) da t a to cosmolog y can be referred 
to dSimon et al.1 (l20C)5[):lYi & Zhand ( l2007h : lsamushia & Ratral ( l2006h : IWei & Zhanj ( l2007h : lLin et al.1 
( l2008l) : lGong et all ( l2008h ) and so on. 



3 CONSTRAINTS ON THE UNDULANT UNIVERSE 

Considering that the contents of the universe are mainly dark matter and dark energy, we just show 
b — ilnio plane contour maps, where b traces dark energy, and fimo traces dark matter respectively. We 
estimate the best fit to the set of parameters by using individual statistics, with 



v2 

A. |data 



E 



[d\iz) - dUz)? 



(8) 



where N is the number of data in each set, and d[{z) is the value at Zi given by above theoretical analysis. 
^0(2;) is individual observational value in each data set at Zi, and o-j is the error due to u ncertainties. 
In ord er to strengthen the constraints, we use MCMC technique to modulate the process dOong et al.l 
(l2007l) ). 

It is worth noting that the value of h may strongly affect the constraint process, because the proba- 
bility distribution function (PDF) of h for each data is different. We can see the PDF for H{z) and SNe 




Fig. 2 The contour maps of 6 — ri,„o for different sets. The left panel is H{z)+BAO, and the 
right panel is SNe la+BAO. Confidence regions are at 68.3%, 95.4% and 99.7% levels from 
inner to outer respectively, for a flat Undulant Universe without a prior of h. 



la data in Fig.5. For H{z) data, we find the PDF of h covers a large range, while it is small for SNe la 
date. Because the distribution of h constrained by SNe la is very narrow, if we fix h far from the best-fit, 
then the constraints will be totally wrong. For example, we test the effect on fixing a "wrong" h by using 
SNe la+BAO+CMB. As can be seen in Fig. 6, with a "wrong" h prior, the constraints of the left panel 
are very different from that of the right panel. This is because if we assume a "wrong" prior for h, the 
would become very large (x^ ^ 426 for /i = 0.65 while ^ 312 for h = 0.7), and then we will get 
a false minimum to lead to a wrong result. However, if we marginalize the likelihood functions over 
h by integrating the probability density P cx 1"^, we can remove the effect of the distribution of ft,, 
and illustrate the reliable distribution of the other parameters. Therefore, we marginalize the likelihood 
functions over h to obtain the best fitting results. Then we get the confidence regions and the h — ilmo 
plane. 

As it is seen from Eq.(l), Loia) is an even function. Correspondingly, the values of h are along the 
axis of ilmo symmetrical distribution. Thus, the constraints on h just in the range from to +oo are 
enough to demonstrate the comparison. Fig. 1 shows the confidence regions determined by the (z) and 
SNe la data alone respectively. We find that the constraints on ri,uo from the Hiz) data are weaker than 
that from the SNe la data. Although we just have 9 data in ll(z) and 307 data in SNe la, there is still 
some consistency between them in the constraints. The joint constraints are shown in Fig. 2, Fig. 3, and 
Fig.4. We find that the constraints are weak if we use the or SNe la data alone, however, after 

inclusion of the other observational data (the BAO or CMB data), the constraints are improved remark- 
ably. Especially, the CMB data improve the constraints more remarkably. For breaking the degeneracy, 
the BAO and CMB data play very important roles in the constraint compared with the H^z) and SNe 
la data alone. In the left contour map of Fig.l, and the two contour maps of Fig. 2, Fig. 3, and Fig.4, 
all contours at the Icr confidence for 6 are open-ended. It shows that 6 takes possible value around 0. 
We have not demonstrated the details of the constraints, because MCMC technique cannot determine 
exact values reliably. These contours correspond to Icr, 2(t, and 3(t confidence from inside to outside. 
Table|2] summarizes 68.3% confidence intervals. From the table, we find that the constraints on ilmo 
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Fig. 3 The contour maps of 6 — fimo for different sets. The left panel is H{z)+CMB, and the 
right panel is SNe la+CMB. Confidence regions are at 68.3%, 95.4% and 99.7% levels from 
inner to outer respectively, for a flat Undulant Universe without a prior of h. 



Table 2 68.3% confidence intervals for each observational data 



Data 


H{z) 


SNe 


H(z)+BAO 


SNe+BAO 


b 


[0.158, 0.43J 
[0,2.0] 


[0.150,0.299] 
[1.13,3.21] 


[0.239,0.315] 
[0, 1.95] 


[0.239, 0.305] 
[0, 2.47] 




Data 


H{z)+CMB 


SNe+CMB 


H{z)+BAO+CMB 


SNe+BAO+CMB 


b 


[0.21,0.30] 
[0, 0.40] 


[0.244,0.318] 
[0, 0.53] 


[0.244, 0.318] 
[0,0.51] 


[0.253, 0.304] 
[0, 0.58] 



from H {z)+BAO are looser than that from SNe la+BAO, and there is little difference between them. 
However, the constraints on b from SNe la+BAO are looser than that from iJ(2;)+BAO. The similar 
conclusions can be derived from comparing H{z)+CMB with SNe la+CMB, and iJ(2)+BAO+CMB 
with SNe la+BAO+CMB. The results of fl y^o on the Undulant U niverse are consistent with the recent 
reported results from WMAP 5-year data ( Kowalski et al.l ( l2008h ). The distribution of b is close to 0, 
i.e. the Undulant Universe approaches the ACDM model. Comparing Fig.4 with Fig. 6, we find that 
the result of the right panel in Fig. 6 with h = 0.7 is similar with that of the right panel in Fig.4 with 
marginalizing h, but has a smaller range for Jlmo. On the other hand, from the result of the left panel 
in Fig. 6 with h = 0.65, we find that the range of b shrinks to 0.6 ^ 0.8 at 68.3% C.L. and the best-fit 
value of ilmo moves to ^ 0.45 around, which are wrong results. 



4 CONCLUSIONS AND DISCUSSIONS 

So far, we've presented some kinds of constraints on cosmological parameters by using the H{z), SNe 
la, BAO and CMB data on the flat Undulant Universe, and we have marginalized the parameter h to get 
the best fitting results. Because if we fix a "wrong" prior for h, any model we take cannot fit the observa- 
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Fig. 4 The contour maps of b— il„^o for different sets. The left panel is iJ (z)+BAO+CMB, and 
the right panel is SNe la+BAO+CMB. Confidence regions are at 68.3%, 95.4% and 99.7% 
levels from inner to outer respectively, for a flat Undulant Universe without a prior of h. 



tional data, and then leads to a false ■ We show confidence regions on the b — Omo plane. The Undulant 
model can trace the distribution of dark energy in undulant case. Therefore, our conclusions are more 
general. As discussed above, we have showed the comparison of different data sets carefully. Comparing 
with the H{z) and SNe la data alone, the BAO and CMB data play very important roles in breaking the 
degeneracy. The constraints on ilnio from the H{z) data are looser than that from the SNe la data, but 
there is little difference between them. For the combined data sets, the constraints are improved remark- 
ably. The CMB data improve the constraints more remarkably than the BAO data. In the joint constraints, 
the results from the H{z) data sets (H{z)+BAO, H{z)+CMB and iJ(z)+BAO+CMB) are almost con- 
sistent with that from the SNe la data sets (SNe la-nBAO, SNe la-nCMB and SNe laH-BAO-nCMB). We 
find all the constraints on b are around 0, but the H{z) data sets constrain b more tightly. Namely, the 
Undulant Universe approaches the ACDM model. However, on ilmo, the constraints of the H{z) data 
sets are weaker than that of the SNe la data s ets. It is probably because the data amount is not sufficient 
and the corresponding errors are very large (ISamushia & Ratral ( |2006|) ). Only are there 9 H{z) data in 
our discussion, but we still get such fine tight constraints. We speculate that the constraints of the H{z) 
data sets will be tighter if we get more H{z) data. Therefore, much more H{z) data have an advantage 
over the SNe la data in the constraint. Fortunately, a large amount of H{z) including data from the AGN 
and Galaxy Evolution Survey (AG ES) and the Atacama Cosm ology Telescope (ACT) are expected to 
be available in the next few years jSamushia & Ratral (l2006l) ). Provided a statistical sample of many 
hundreds of galaxies, we could determine the value of H{z) data to a percent accuracy. Paying attention 
to the assumption that SNe la is "standard candle", it might make the luminosity distance of SNe 
la sample biased. In addition, there exists an integration of the inverse of H{z) in d^, and some uncer- 
tainties in the minimizing statistics might also arise from this integration. On the contrary, H{z) do 
not suffer from this uncertainty resulting from integration. Therefore, the H{z) data should constrain 
parameters more strongly than the SNe la data. We also notice that, the redshift range of H{z) is from 
0.09 to 1.75 and the redshift range of SNe la is from about 0.01 to 1.75. Thus, H{z) and SNe la are 
distributed in the same range of redshift roughly. Therefore, the H{z) data seem to provide us with 
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Fig. 5 The probability distribution function of h for H{z) and SNe la data, respectively, for 
a flat Undulant Universe. 

an independent, very simple, and very powerful probe of fundamental cosmological parameters. With 
a large amount of the H{z) data in the future, we probably can constrain cosmological parameters by 
using the H{z) data instead of the SNe la data. 
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